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66 ■ ABSTRACT 

o ■ 
o : 

CN . Mergers between stellar-mass black holes will be key sources of gravitational 

p • radiation for ground-based detectors. However, the rates of these events are 

■ highly uncertain, given that such systems are invisible. One formation scenario 

. involves mergers in field binaries, where our lack of complete understanding of 

^ . common envelopes and the distribution of supernova kicks has led to rate es- 

i"^. timates that range over a factor of several hundred. A different, and highly 

Q^. promising, channel involves multiple encounters of binaries in globular clusters 

Q . or young star clusters. However, we currently lack solid evidence for black holes 

. in almost all such clusters, and their low escape speeds raise the possibility that 

ci . most are ejected because of supernova recoil. Here we propose that a robust 

environment for mergers could be the nuclear star clusters found in the centers 

^ . of small galaxies. These clusters have millions of stars, black hole relaxation 

^ . times well under a Hubble time, and escape speeds that are several times those 

, of globulars, hence they retain most of their black holes. We present simulations 

^ ! of the three-body dynamics of black holes in this environment and estimate that, 

^ ! if most nuclear star clusters do not have supermassive black holes that interfere 

OO ! with the mergers, at least several tens of events per year will be detectable with 

'. Advanced LIGO. 
> ■ 

^ ■ Subject headings: black hole physics - galaxies: nuclei - gravitational waves — 
' 

' relativity 



1. Introduction 



Groun d-based gravitationa l wave detectors have now achieved their initial sensitivity 



goals (e.g., lAbbott et al.ll2007l ). In the next few years, these sensitivities are expected to 
improve by a factor of ~ 10, which will increase the searchable volume by a factor of ~ 10^ 
and will lead to many detections per year. 



One of the most intriguing possible sources for such detectors is the coalescence of a 
double stellar-mass black hole binary. Such binaries are inherently invisible, meaning that 
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we have no direct observational handle on how common they are or their masses, spin 
magnitudes, or orientations. Comparison of the observed waveforms (or of waveforms from 
merging supermassive black holes) with predictions based on post-Newtonian analysis and 
numerical relativity will be the most direct possible test of the predictions of strong-gravity 
general relativity. 

The electromagnetic non-detection of these sources makes rate estimates highly chal- 
lenging, because our only observational handles on BH-BH binaries come from possible 
progenitors. For example, a common scenario involves the effectively isolated evolution of 
a field binary containi ng two massive star s into a binary with two black holes that will 



eventually merge (e.g., iLipunov et al.l 119971 : iBelczynski &: Buliklll999l ). There are profound 



uncertainties involved in calculations of these rates due to e.g., the lack of knowledge of the 
details of the common envelope phase in these systems and the absence of guides to the 
distribution of supernova kicks delivered to black holes. As a recent indication of the range 
of estimated rates, note that the Advanced LIGO d etection rate of BH-BH coalescences is 
estimated to be anywhere between ~ 1 — 500 yr~^ by IBelczynski et al.l (l2007l ). depending on 
how common envelopes are modeled. 

Another promising location for BH-BH mergers is globular clusters or super star clusters, 
where stellar number densities are high enough to cause multiple encounters and hardening of 



binaries. Even though binaries are kicked out before they nierge ( Kulkarni. Hut, fc McMillan 



1993 



2000 



O'Leary et al. 



Sigurdsson &: Hernauistlll993l : ISigurdsson &: Phinneylll993 



1995 



Portegies Zwart &: McMillan 



200 



these clusters can still serve as breeding grounds for gravitational 



O'Leary. O'Shaughnessy. fc Rasiol (l2007l ) estimate a rate of 0.5 yr 



wave sources. Indeed, 

for initial LIGO and 500 yr~^ for Advanced LIGO via this channel. There is, however, little 
direct evidence for black holes in most globulars (albeit they could be difficult to see). In ad- 
dition, at least one black hole in a low-mass X-ray binary apparen tly received a ^ 100 km s"^ 
kick from its supernova (GRO J1655-40; see lMirabel et al.l 120021). This is double the escape 
speed from the centers of even fairly rich globulars JWebbinklll985[ ). leading to uncertainties 
about their initial black hole population and current merger rates. 

Here we propose that mergers occ ur frequently in the nuclear star clusters that may be in 
the c enters of many low- mass galaxies ( Boker et al.ll2002 : Ferrarese et al.|[2006 : Wehner fc Harris 
20061 : note that some of these are based on small deviations from smooth surface bright- 



ness profiles and are thus still under discussion). It has recently been recognized that 
in these galaxies, which may not have supermassive bl ack holes (for a sta tus report on 
ongoing searches for low-mass central black holes, see iGreene fc Hoi 120071 ). the nuclear 
clusters have masses that are correlated with the surrounding velocity dispersion a as 
M 10^ M0(cr/54 km s~^)^'^ ( Ferrarese et al.ll2006 ). When the velocity dispersion is in 
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the range of ~ 30 — 60 km s~^, the half-mass relaxation time is small enough that black 
holes (which have ~ 20x the average stellar mass) can sink to the center in much less than 
a Hubble time. In addition, although systems wit h equal-mass objects requ ire roughly 15 
half-mass relaxation times to undergo core collapse (IBinney fc Tremainelll987l ). studies show 
that systems with a wide r ange of stellar masses experierice core collapse within ~ 0.2 x the 
half-mass relaxation time (jPortegies Zwart fc McMillanI 120021 : iGiirkan et al.l |2004| ) . There- 
fore, clusters with masses less than ~ few x 10^ Mq will have collapsed by now and hence 
increased the escape speed from the center, allowing retention of most of their black holes. 

As we show in this paper, nuclear star clusters are therefore excellent candidates for 
stellar-mass black hole binary mergers because they keep their black holes while also evolving 
rapidly enough that the holes can sink to a region of high density. If tens of percent of the 
black holes in eligible galaxies undergo such mergers, the resulting rate for Advanced LIGO is 
at least several tens per year. In § 2 we quantify these statements and results more precisely 
and discuss our numerical three-body method. We give our conclusions in § 3. 



Method and Results 



2.1. Characteristic times and initial setup 



Our approach is similar to that of lO'Leary et al.l (120061 ) . who focus on globular clusters 
with velocity dispersions a < 20 km s~^. Here, however, we concentrate on the more 
m assive and tightly boun d nuclear star clusters. Our departure point is the relation found 
by iFerrarese et al.l (120061 ) between the masses and velocity dispersions of such clusters: 



10' 



6.91±0.11 



((7/54 km s ^) 



1\4.27±0.61 



Assuming that there is no massive central black hole for these low velocity dispersions, the 



half-mass relaxation time for the system is (see iBinney fc Tremaind Il987l ) trix ~ g^^f Across 



^nuc/0.5 Mq is the number of stars in the system (assuming an average mass 



where 

of 0.5 Mq) and tcmss ~ -R/cr is the crossing time. Here R 
cluster. Putting this together gives 



GMnuc/o" is the radius of the 



trix ~ 1-3 X 10^° yr ((t/54 km s~^) 



n 5.54±0.61 



(2) 



The r elaxation time scales inversely with the mass of an individual star (IBinney &: Tremaine 
19871 ). so a 10 Mq black hole will settle in roughly 1/20 of this time. Also note that large 
N-body simulations with broad mass functions evolve to core collapse within ro ughly 0.2 
half-mass relaxation times (jPortegies Zwart fc McMillarul2002l : iGiirkan et al.ll2004l ). hence in 
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the current universe clusters with velocity dispersions a < 60 km s ^ will have had their 
central potentials deepened significantly. 

The amount of deepening of the potential, and thus the escape speed from the center of 
the cluster, depends on uncertain details such as the initial radial dependence of the density 
and the binary fraction. Given that the timescale for segregation of the black holes in the 
center is much less than a Hubble time, we will assume that the escape spe ed is roughly 5x 
the velocity dispersion, as is the case for relatively rich globular clusters (jWebbinkI Il985l ) . 
This may well be somewhat conservative, because the higher velocity dispersion here than in 
globulars suggests that a larger fraction of binaries will be destroyed in nuclear star clusters. 
This could lead to less efficient central energy production and hence deeper core collapse 
than is typical in globulars. 

With this setup, our task is to follow the interactions of black holes in the central regions 
of nuclear star clusters, where we will scale by stellar number densities of n ~ 10^ pc~^ 
because of density enhancements caused by relaxation and mass segregation. Some black 
holes will begin their lives in binaries, but to be conservative we will assume that they start 
as single objects and have to exchange into binaries that contain main sequence stars or 
other objects. All binaries in the cluster will be hard, i.e., will have internal energies greater 
than the average kin etic energy of a field star, because otherwise they will be softened and 



ionized quickly (e.g., iBinney fc Tremaindll987l ). If, for example, we consider binaries of two 



1 Mq stars in a system with ex = 50 km s~^, this means that the semimajor axis has to be 
less than Omax ~ 1 AU. Studies of main sequence binaries in globular clusters, which have 
cr ~ 10 km s~^, suggest that after billions of years roughly 5-20% of them survive, with the 
rest falling victim to ionization or collisions (jivanova et al.l 120051 ). The binary fraction will 
be lower in nuclear star clusters due to their enhanced velocity dispersion, but since when 
binaries are born they appear to have a con s tant distribution across the log of the semimajor 
axis from ~ 10^^ — 10^ AU (e.g., lAbtlll983l : iDuquennoy fc Mayoiill99ll ) the reduction is not 
necessarily by a large factor. We will conservatively scale by a binary fraction /bin = 0.01. 

If a black hole with mass Mbh gets within a couple of semimajor axes of a main se- 
quence binary, the binary will tidally separate and the BH will acquire a companion. The 



timescale on which this happens is t 



bin = {nT.a)-\ where S = vrr^ (1 + 2GMtot/(aVp)) is 
the interaction cross section for pericenter distances < Tp when gravitational focusing is in- 
cluded. Here Mtot is the mass of the black hole plus the mass of the binary. If we assume 
that Mbh = 10 Mq and it interacts with a binary with two 1 Mq members and an a = 1 AU 
semimajor axis, then the typical timescale on which a three-body interaction and capture of 
one of the stars occurs is 



ts-bod = {nY.a)-^ ^ 1.2 x 10^ yr(n/10^ pc-3)-i(/bin/0.01)-^(cT/50 km s)(a/l AU)"^ . (3) 
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This is small enough compared to a Hubble time that we start our simulations by assuming 
that each black hole has exchanged into a hard binary, and follow its evolution from there. 

Another important question is whether, after a three-body interaction, a black hole 
binary will shed the kinetic energy of its center of mass via dynamical friction and sink to 
the center of the cluster before another three-body encounter. If not, the kick speeds will 
add in a random walk, thus increasing the ejection fraction. 

To compute this we note that the local relaxation time of a binary is 

~ In A G^m)Munn ^ ^ 



( ISpitzerl 119871 ) where a is the local velocity dispersion. In A ~ 10 is the Coulomb logarithm, 
(m) is the average mass of interloping stars, n is their number density, and Mbin is the mass 
of the binary. The timescale for a three-body interaction is ts-bod = (nSa)"^ as above. For 
a gravitationally focused binary, which is of greatest interest because only these could in 
principle produce three-body recoil sufficient to eject binaries or singles, Vp < GMbm/o"^. 



Therefore, S ^ 27rrpG'Afbin/o'^ and 



f-S-bod 



27mrpGMua 
If we let Vp = gGMbin/o"^, with q < 1, then 



(5) 



f3-bod 



so that 



(6) 



t.l./t3-bod - • (7) 

In the center of a cluster, where mass segregation is likely to have flattened the mass distri- 
bution, we flnd that this quantity is typically less than unity (and it decreases as the binary 
hardens), meaning that after a three-body encounter a binary has an opportunity to share 
its excess kinetic energy via two-body encounters and thus settle back to the center of the 
cluster. We therefore treat the encounters separately rather than adding the kick speeds in 
a random walk. 

In a given encounter, suppose that a binary of total mass Mbm = Mi + M2, a reduced 
mass fi = MiM2/Mbin, and a semimajor axis Oinit interacts with an interloper of mass mint, 
and that the kinetic energy of the interloper at inflnity is much less than the binding energy 
of the binary (i.e., this is a very hard interaction). If after the interaction the semimajor axis 
is ttfin < ^init, then energy and momentum conservation mean that the recoil speed of the 
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binary is given by f^jn = ^/^ Mb"+'m- t (V'^fin ~ l/'^init), and the recoil speed of the interloper 
is Vint = (Afbin/'^int)'i^bin- For example, suppose that Mi = M2 = 10 Mq, Mint = 1 M©, 
Oinit = 0.1 AU, and afin = 0.09 AU. The binary then recoils at fbin = 15 km s~^ and stays in 
the cluster, whereas the interloper recoils at Vi^t = 300 km s~^ and is ejected. 



2.2. Results 



The central regions of the clusters undergo significant mass segregation, and thus the 
mass function will be a t least flatt ened, and possibly inverted. This has, for example, been 
observed for globulars (jSosinI 119971 ). To include this effect, when we consider the mass of a 
black hole, its companion, or the interloping third object in a binary-single encounter, we go 
through two steps. First we select a zero age main sequence (ZAMS) mass between 0.2 Mq 
and 100 Mq using a simple power law distribution dN/dM oc M~". We allow a to range 
anywhere from 2.35 (the unmodified Salpeter distribution) to —1.0, where smaller values 
indicate the effects of mass segregation. Second, we evolve the ZAMS mass to a current 
mass. Our mapping is that for Mzams < 1 Mq, the star is still on the main sequence and 
retains its original mass; for 1 Mq < Mzams < 8 Mq the star has evolved to a white dwarf, 
with mass MwD = 0.6 Mq + 0.4 Mq(Mzams/M0 - 0.6)^/3. 8 Mq < Mzams < 25 M© 
the star has evolved to a neutron star, with mass M^s = 1-5 Mq + 0.5 M© (Mzams ~ 
8 M0)/17 Mq] and for Mzams > 25 Mq the star has evolved to a black hole with mass 
Mbh = 3 M0 + 17 M0 (Mzams - 25 Mq) /75 Mq. Therefore, we assume that black hole 
masses range from 3 Mq to 20 Mq. 

These prescriptions are overly simplified in many ways. We therefore explore different 
mass function slopes, main sequence cutoffs, and so on, and find that our general pic ture 
is robust against specific assumptions. Note that, consistent with lO'Leary et al.l (120061 ) . we 
find that there is a strong tendency for the merged black holes to be biased towards high 
masses. Therefore, if black holes with masses > 20 Mq are common, these will dominate 
the merger rates. This is important for data analysis strategies, because the low-frequency 
cutoff of ground-based gravitational wave detectors implies that higher-mass black holes will 
have proportionally more of their signal in the late inspiral, merger, and ringdown. 

The three-body interactions themselves are assumed to be Newtonian interactions be- 
tween point masses and are computed using the hierarchical N-body code HNBody (K. Ranch 
and D. Hamilton, in preparation), using t he driver lABL developed by Kayhan Giiltekin (see 
Giiltekin. Miller, fc Hamiltonll2004l . |2006| for a detailed description). These codes use a num- 
ber of high-accuracy techniques to follow th e evolution o f gravitating point masses. Between 
interactions, we use the Peters equations (jPeterd Il964j ) to follow the gradual inspiral and 
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circularization of the binary via emission of gravitational radiation. This is neghgible except 
near the end of any given evolution. 

We begin by selecting the mass of the black hole and of its companion (which does not 
need to be a black hole) from the evolved mass function. We also begin with a semimajor 
axis that is 1/4 of the value needed to ensure that the binary is hard. We do this because 
soft binaries are likely to be ionized and thus become single stars rather than merge. We also 
select an eccentricity from a thermal distribution P{e)de = 2ede. We then allow the binary 
to interact with single field stars drawn from the evolved mass function, one at a time, until 
either (1) the binary merges due to gravitational radiation, (2) the binary is split apart and 
thus ionized (this is exceedingly rare given our initial conditions), or (3) the binary is ejected 
from the cluster. The entire set of interactions until merger typically takes millions to tens 
of millions of years, and only rarely over a hundred million years, so it finishes in much 
less than a Hubble time. In the course of these interactions there are typically a number of 
exchanges, which usually swap in more massive for less massive members of t he binary. This 



i s the cause of the bias towards high-mass mergers that was also found by lO'Leary et al. 



(I2OO6I ). As shown in Table 1, for a < 1 most black holes acquire a black hole companion in 



the process of exchanges, and for a < 0.5 virtually all do. 

The results in Table 1 are focused on different mass function slopes and escape speeds. 
As expected, we find that for Vesc > 150 km s~^ the overwhelming majority of black hole 
binaries merge in the nuclear star cluster rather than being ejected (see Figure 1). This is 
the difference from lower-cr globular clusters, where the mergers happen outside the cluster. 
Note also that in addition to few binaries being ejected, there are typically only 1-2 single 
black holes ejected per merger, suggesting that > 50% of holes will merge. In contrast, at the 
50 km s~^ escape speed typical of globulars, > 20 single black holes are ejected per merger, 
suggesting an efficiency of < 10%. For well-segregated clusters (with a < 0), the average 
mass of black holes that merge, binary ejection fraction and number of singles ejected, and 
number of black holes that merge with each other instead of other objects are all insensitive 
to the particular mass function slope. For less segregated clusters with a > 0, the retention 
fraction of black holes rises rapidly to unity because most of the objects that interact with 
the holes are less massive stars. In such cases there might be a channel by which the mass 
of the holes increases via accretion of stars, but we expect a > to be rare for nuclear 
star clusters because of the shortness of the segregation times of black holes. Overall, there 
appears to be a wide range of realistic parameters in which fewer than 10% of binary black 
holes are ejected before merging. 
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Discussion and Conclusions 



We have shown that nuclear star clusters with velocity dispersions around o" ~ 30 — 
60 km s~^ are promising breeding grounds for stellar-mass black hole mergers. At signif- 
icantly lower velocity dispersions, as found in globulars, the escape speed is low enough 
that the binaries are ejected before they merge. Significantly high er velocity dispersions 



appear correlated with th e appearance of supermassive black holes (iGebhardt et al.l 12000 



Ferrarese fc MerrittI l2000l ) . In such an environment there might also be interesting rates of 
black hole mergers, but the increasing velocity dispersion closer to the central object means 
that binary fractions are lower and softening, ionization, or tidal separa tion by the super 



massive black hole itself are strong possibilities for stellar-mass binaries (IMiller et al.l 12005 
Lauburg & Miller, in preparation). 

To estimate the rate of detections with Advanced LIGO, we note that velocity disper- 
sions in the a ~ 30 — 60 km s~ ^ range correspond to roughly a factor of ~ 10 in g alaxy 
luminosity (IFerrarese et al.ll2006l ). Galaxy surveys suggest (e.g., iBlanton et al.ll2003l ) that 
for dim galaxies the luminosity function scales as roughly dN/dL = 0*(L/L*)~", where 
0* = 1.5 X 10~'^h^ Mpc~^ ^ 5 X 10"^ Mpc"^ for h = 0.71, and a ^ —1. This implies 
that there are nearly equal numbers of galaxies in equal logarithmic bins of luminosity. A 
factor of 10 in luminosity is roughly e^, so the number density of relevant galaxies is ap- 
proximately 10^^ Mpc~'^. To get the rate per galaxy, we note that typical initial mass 
functions and estimates of the mass needed to evolve into a black hole combine to suggest 
that f or a cluster of mass M„uc., approximately 3 x lO~'^(Mnuc/^0) stars evolve into black 
holes (lO'Leary. O'Shaughnessy. &: Rasid 120071 ). That implies a few x 10^ black holes per 
nuclear star cluster. If a few tens of percent of these merge in a Hubble time, and if the 
rate is slightly lower now because many of the original black holes have already merged (see 
O'Leary et al. 2006), that suggests a merger rate of > 0.1 x few x 10^/(10^° yr) per galaxy, or 
fewxlO"^ Mpc~^ JT^^. At the ~ 2 Gpc distance at w hich Advanced LIGO is expected to be 



able to see mergers of two 10 Mq black holes (see, e.g.. lMandelll2007l ). the available volume is 
3 X 10^° Mpc^, for a rate of ^ 100 per year. Roughly 5 0-80% of galaxies in the eligible lumi- 
nosity range appear to have nuclear star clusters (see IFerrarese et al.ll2006l for a summary). 
If the majority of the clusters do not have a supermassive black hole, this suggests a final 
rate of at least several tens per year for Advanced LIGO. This could be augmented somewhat 
by small galaxies that originally had superm assive black holes, but had them ejected after a 
merg er and then reformed a central cluster (IVolonteril l2007l : IVolonteri. Haardt. &: Giiltekin 

mm- 



For nearby {z < 0.1) events of this type it might be possible to identify the host 
galaxy. However, for more typical z ~ 0.5 =^ d ^ 2 Gpc events the number of can- 
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didates is too large: even assuming angular localization of AQ = (1°)^ and a distance 
accuracy of Ad/d = 1%, the number of galaxies in the right luminosity range is ~ 
47r/3(2000 Mpc)3(Al]/47r)(Ad/d)(0.01 Mpc"^) ^ 80. Therefore, barring some unforseen 
electromagnetic counterpart, the host will usually not be obvious. 

We anticipate that tens per year is a somewhat conservative number, because our sim- 
ulations suggest that more like 50% of black holes will be retained, even as single objects, 
and because (unlike in a globular cluster) the central regions of galaxies are not devoid of 
gas, hence more black holes could form in the vicinity of the cluster and fall in. In addi- 
tion, if stellar-mass black holes with masses beyond 20 Mq are common, this also increases 
the detection radius and hence the rate. Even for total masses ~ 30 Mq and at redshifts 
z ~ 0.5, the observer frame gravitational wave frequency at the innermost stable circular 
orbit is /isco ~ 4400 Hz/[30(1 + z)] ~ 100 Hz. This is close enough to the range where 
frequency sensitivity declines that detection of many of these events will rely strongly on the 
signal obtained from the last few orbits plus merger and ringdown. In much of this range, 
numerical relativity is essential. 

As a final point, we note that for the same reason that nuclear star clusters are fa- 
vorable environments for retention and mergers of stellar-mass black holes, they could 
also be good birthplaces for more massive black holes. This could be prevented, even 
for the relatively high escape speeds discussed here, if recoil from gravitational radiation 
during the coalescence exceeds ~ 200 km s~^. The key uncertainty here is the spin mag- 
nitudes of the holes at birth. Numerous simulations demonstrate that high spins with sig- 
nificant projections in the bi nary orbital plane ca n produce kicks of up to several thou- 
sand kilometers per second (j Gonzalez et al.l 120071 ). If there is significant processing of 
gas through accretion disks the spins are aligned in a way that reduces the kick to below 
200 km s"^ flBogdanovic. Revnolds. fc Millei]l2007f ). but stellar -mass black holes cannot pick 
up enough mass from the interstellar medium for this to be effective. For example, the Bondi- 
Hoyle accretion rate is MBondi ~ 10"^^ Mq yr-i((T/50 km s~^)-^{n/100 cm-3)(M/10 Mq)^, 
meaning that to accrete the ~ 1% of the black hole mass needed to realign the spin 



(IBogdanovic. Reynolds, fc Millerl 120071 ) would require at least a trillion years. Current es- 



timates of stellar-mass black hole spins suggest a/M > 0.5 in many cases (IShafee et al. 



20061 : iMcClintock et ahlbood : lMillerll2007l : lliu et al.l 120081 ). If the spins are isotropically ori- 
ented and uniformly distributed in the range < a/M < 1, and the mas s ratios are in the 



/mh| g ~ 0.6 — 0. 8 rang e typical in our simulations, then use of the ICampanelli et al. 



(120071 ) or iBaker et al.l (120081 ) kick formulae imply that roughly 84% of the recoils exceed 
200 km s~^ and 78% exceed 250 km s~^. This suggests that r nultiple mergers are rare unless 



there is initially an extra-massive black hole as a seed (e.g., iHoUey-Bockelmann et al.l 12008 



for a discussion of the effects of gravitational wave recoil), but further study is important. 



-lo- 



in conclusion, we show that the compact nuclear star clusters found in the centers of 
many small galaxies are ideal places to foster mergers between stellar-mass black holes. It is 
not clear whether multiple rounds of mergers can lead to runaway, but this is a new potential 
source for ground-based detectors such as Advanced LIGO, where numerical relativity will 
play an especially important role. 

We thank Karl Gebhardt, Kayhan Giiltekin, Kelly HoUey-Bockelmann, and Jeff Mc- 
Clintock for their many useful suggestions, which helped to clarify this paper. This work 
was supported in part by NASA ATFP grant NNX08AH29G. 
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Table 1. Simulations of nuclear star clusters" 



,max 

,merge e,eject/ 



50 


IMo 





11.7 


62.5 


IMq 





11.7 


75 


IMq 





11.7 


87.5 


IMq 





11.7 


20 


IMq 





11.7 


112.5 


IMq 





11.7 


125 


IMq 





11.7 


137.5 


IMq 





11.7 


150 


IMq 





11.7 


162.5 


IMq 





11.7 


175 


IMq 





11.7 


187.5 


IMq 





11.7 


200 


IMq 





11.7 


212.5 


IMq 





11.7 


225 


IMq 





11.7 


237.5 


IMq 





11.7 


250 


IMq 





11.7 


200 


IMq 


-1.0 


13.4 


200 


IMq 


-0.5 


12.6 


200 


IMq 


0.5 


10.7 


200 


IMq 


1.0 


9.7 


200 


IMq 


1.5 


8.8 


200 


IMq 


2.0 


7.5 


200 


IMq 


2.35 


7.4 


200 


3Mq 


-1.0 


13.4 


200 


3Mq 


-0.5 


12.6 


200 


3Mq 





11.7 


200 


3Mq 


0.5 


10.7 


200 


3Mq 


1.0 


9.7 


200 


3Mq 


1.5 


8.8 


200 


3Mq 


2.0 


7.5 



0.25 


0.0 


31.2 


24.8 


0.33 


0.0 


31.6 


15.3 


0.42 


0.0 


30.9 


11.5 


0.52 


0.0 


31.9 


7.9 


0.63 


0.02 


30.0 


6.2 


0.68 


0.0 


31.4 


4.7 


0.72 


0.02 


31.8 


4.3 


0.76 


0.01 


32.0 


3.0 


0.80 


0.03 


32.3 


2.8 


0.93 


0.03 


31.3 


2.0 


0.89 


0.02 


31.9 


2.0 


0.90 


0.01 


31.3 


2.1 


0.94 


0.08 


31.1 


1.3 


0.89 


0.05 


30.5 


1.0 


0.98 


0.06 


31.0 


1.2 


0.94 


0.06 


30.1 


1.0 


0.96 


0.06 


30.0 


0.71 


0.94 





32.4 


1.3 


0.95 


0.01 


32.2 


1.5 


0.94 


0.1 


28.3 


0.91 


0.98 


0.41 


27.3 


0.43 


0.99 


0.79 


23.0 


0.04 


1.00 


0.99 







1.00 


1.00 







0.85 


0.03 


33.3 


1.5 


0.94 


0.01 


31.9 


1.3 


0.95 


0.05 


30.4 


1.5 


0.94 


0.11 


29.2 


1.0 


0.99 


0.48 


25.3 


0.38 


0.99 


0.85 


24.7 


0.04 


1.00 


1.00 
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Table 1- 


—Continued 






K,, (km l/,ns.nK.x^- n'' ( 


A^ii)(-^/. 






\f . . V' 

^' siii,i^i(\t\i('<.*i7 


200 3Mq 2.35 


7.4 


1.00 1.00 — 







*A11 runs had 100 realizations. 



''Escape speed from cluster. 
'^Maximum mass of main sequence star. 

"^Number distribution of stars on zero age main sequence: dN/dM oc M~". 

^Average mass of all black holes given a and our evolutionary assumptions. 

^Fraction of runs in which holes merged rather than being ejected. 

^Fraction of runs in which holes merged with something other than another black hole. 

'^Average mass of double BH binaries that merged. 

'Average number of single black holes ejected per binary that merged. 
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Fig. 1. — Fraction of binaries retained in the nuclear star cluster (solid line) and average 
number of black holes ejected per black hole merger (dotted line) as a function of the cluster 
escape speed. Here the zero age main sequence distribution of masses is dN/dM oc M'^, to 
account for mass segregation in the cluster center, where most interactions occur. We also 
assume a maximum black hole mass of 20 Mq and a maximum main sequence mass of 1 M©, 
but most results are robust against variations of these quantities. All runs are done with 100 
reahzations, which explains the lack of perfect smoothness. We see, as expected, that the 
retention fraction increases rapidly with escape speed, so that for nuclear star clusters most 
binaries stay in the cluster until merger. We also see that at Ksc ~ 200 km s~^ and above, 
most black hole singles also stay in the cluster. This suggests a high merger efficiency. 



